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Abstract 

We study high-prHiggs boson production at hadron colliders to order 0(a s Gp) in 
hadron collisions. In particular, we investigate the process g + q/q — > q/q + H, where 
q = u,d,c,s, or b, for the LHC (a y/s =14 TeV, proton-proton collider). Our results 
are compared to the ©(o^G^) calculation. The associated production of a high-py 
Higgs boson with a 6-quark or anti-quark is comparable to the 0(a^Gj?) calculation 
because of the large top quark mass and the additional contribution of electroweak 
gauge and Goldstone bosons. The associated production of light quarks, however, is not 
significant. We also comment on new physics effects in the framework of the electroweak 
chiral Lagrangian. 
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1 Introduction 



With the discovery of the top quark the only remaining element of the Standard 
Model (SM) particle spectrum is the Higgs boson. Experimentally, there are only lower 
bounds on Mjj. LEP-I has placed the limit Mjj > 64.5 GeV 0. Theoretically, there are 
upper bounds in the SM from unitarity and triviality arguments |Q . One goal of the future 
High Energy Physics experimental program is to discover the Higgs boson and verify its 
properties or determine the alternative mechanism of electroweak symmetry breaking. The 
search for the Higgs boson at LEP-II is strictly limited by the available center of mass 
energy and luminosity, so that only Mg < 90-95 GeV can be probed for an energy of 
190 GeV and 500 pb _1 of data 0. The reach of a high-luminosity Tevatron collider is 
better, but becomes challenging above Mjj = 110 GeV |j. The LHC, on the other hand, 
is hoped to have enough energy, luminosity, and instrumentation to decisively probe the 
energy scale associated with electroweak symmetry breaking. This task is not as straight- 
forward as it may seem. One possible alternative to the SM is the Minimal Supersymmetric 
Standard Model (MSSM) with a constrained multi-dimensional parameter space [||, 0, g. 
The constrained MSSM models predict that the couplings of the lightest Higgs boson to SM 
particles is very SM-like, i.e. sin 2 (/?— a) ~ 1, and its mass should be less than about 140 GeV 
|]]. In this case, to deduce supersymmetry, one must observe a superpartner directly or 
discern its presence in quantum corrections. Another alternative, a model with a strongly 
interacting scalar sector [[!(]] , predicts a greatly enhanced Higgs boson width even for a Higgs 
mass of a few hundred GeV, so that the Higgs boson signal can be hidden by backgrounds. 
Regardless of the scenario, a full verification of the properties of the Higgs boson requires a 
deep theoretical understanding of its properties. In this investigation, we concentrate on the 
high-pTproduction of Higgs bosons, which is sensitive to loop corrections. The 0{a ? i.Gp) 



contribution to high-pj-Higgs boson production was calculated previously [11|, where Gp = 
(VZv 2 )" 1 and the vacuum expectation value v = 246 GeV. Here, we extend that calculation 
to include the 0(a s G F ) contributions from electroweak gauge bosons, Goldstone bosons, 
and quarks. In particular, since the top mass is large, we expect to see an enhancement in 
the associated production of a Higgs boson with a 6-quark or anti-quark in some kinematic 
region. We also expect this channel to be sensitive to the coupling of the electroweak 
gauge bosons and Goldstone bosons to the Higgs boson, since it does not vanish in the 
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limit that the U(l)y and SU(2)l gauge couplings vanish.f] We study this sensitivity in 
the framework of the electroweak chiral Lagrangian, which allows us to construct the most 
general effective Lagrangian that is consistent with SU{2)l x U(1)y — > U(l) em symmetry 
breaking. We show that with new physics, the 0(a s Gp) contribution can be comparable 
to the 0{oft s Gp) contribution for high-p^Higgs boson production. Another source of high- 
p^Higgs bosons is the 0(a s Gp) tree level process. We show that this is small by examining 
the processes q + q — > b + b + H and g + g — > b + b + H. We also argue that any interference 
between this order amplitude and one of higher order is suppressed because the bottom 
quark mass m& is much less than the electroweak symmetry breaking scale v. 

2 High— ^production of the Higgs boson to 0(a s Gp) 

To 0(a s G F ), the Higgs boson is produced at high-p^from quark-gluon, antiquark- 
gluon, and quark-antiquark initial states. Since they are the most interesting, we will con- 
centrate on the first two processes for the purpose of this discussion. The quark-antiquark 
annihilation process is typically an order of magnitude smaller at the LHC for Mh < 400 
GeV. We chose to perform the calculation in the helicity formalism, since the amount of 
algebra is reduced significantly. Furthermore, we used the Feynman rules in the 't Hooft- 
Feynman gauge, since the electroweak gauge bosons and their associated Goldstone bosons 
have the same mass and, hence, loop integrals involving gauge bosons and Goldstone bosons 
have the same denominators. This choice is also advantageous for investigating the elec- 
troweak chiral Lagrangian. 

We consider the process g(p g ) + q(pj) — ► q(Pi) + H(ph) to 0(a s Gp), where the p g and pj 
are the four-momenta of the incoming particles and Pi and pu are the four-momenta of the 
outgoing particles. The quark-antiquark initial state can be generated by the substitution 
Pi — > —pi,p g — * —p g and a reevaluation of the color factor. Contributions to the loop 
integral come from internal lines involving the weak isospin quark partner of q (we use 
the simplification V u d = V cs = Vtb = 1), gauge bosons, and Goldstone bosons. Some 
representative Feynman diagrams are illustrated in Figure 1. The amplitude for this process 

3 The 0(o!sGf) contribution does not depend on the electroweak gauge couplings, but is only sensitive 
to the coupling of the top quark to the Higgs boson. 
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Figure 1: Some representative Feynman diagrams for the C(a s G ! |') contributions to the 
process g + 6/6 — ► 6/6 + H . 

can be written as: 

M\ lXjXg = i9a 1 B()*,Pi)[J : 'yp + grfB]v(\j,P3)f% g , (1) 

where Aj and Xj are the fermion helicity indices, X g is the gluon polarization, 7^ are 4x4 
gamma matrices, ^ = a M 7^ for 4-vector a^, and equations of motion have been applied to 
the on-shell, 4-component spinors. 

In addition, we have taken the limit m q — * for all q except t. The complex scalar T and 
the complex vector are form factors resulting from the integration of the loop momentum, 
and their explicit expression are given in the Appendices I and II. The helicity amplitude 
can be re-written in terms of 2-component Weyl spinors using the bra-ket notation^ 



U-{\ = -1/2, p a ) = w+|Pa-}, 

V_(\ = +1/2, Po ) = -<4|pa-), 

u+(A = +1/2, p a ) = c4|p a +), 

T7+(A = -l/2,p ) = -<|pa+>, 

where u)°L = y/2E a for massless fermions with energy E a , and (p ± | = (|p±))'. These 

are the only components for massless fermions. For p^ = (Ejpsec^jpsgs^jpcg), where 

s$ and C0 are shorthand for sin-0 and cosip, |p+) = (cos #/2, e l< ^ sin #/2) T , and |p— ) = 
4 The 2-component Weyl spinors are denned by the relation u± = |(1 ± 75)11, etc. 
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(— e % ^ sin 0/2, cos #/2) T []. Also, the gluon polarization 4-vectors for left-handed (L) and 
right-handed (i?) helicities can be written as: 



e -i<f> 

e^ L) = — -=-[0, is,/, + c^cg, —ic<f, + s^ce, -so], 

e i4> 

where and 6 are spherical coordinates of the gluon momentum. In the helicity basis, there 
are four non-vanishing amplitudes (as m q — > 0). The parton-level cross section is: 

da = j±C gb £ (\M^ x f + \M ++x f)dR 2 , 

M—\ g = ig 8 u+v+fa - WFl+p + Gn$H + ]\Pj-)£x g i 
M ++Xg = igstj+fJ+fc + |[^7- M + G^H-]\Pj+)e^ g , 

where the flux factor F = 2s for s = (p g + Pj) 2 ', the spin average factor 5 = 2x2; the 
color factor is C gq = 4/(3 x 8) for g + q — ► q + H, and C gg = 4/(3 x 3) for g + q — > # + #; 
the gluon polarization is specified by A 9 ; is the two-body phase space; 7± are the 
2x2 matrices (1, T<?i)j[] and </l± = a^7±. In the above result, M. ++ \ g , which only contains 
contributions from Z°-bosons, are small for two reasons. First, since we are only interested 
in initial and final states without t-quarks, the internal quark is always light (because of the 
neutral current) and has a tiny coupling to the Goldstone boson. Secondly, the left- and 
right-handed couplings of the Z°-boson, which are smaller than the purely left-handed 
coupling of the VF^-bosons, appear in the squared matrix element to the fourth-power. 
For all practical purposes, then, the Z° contributions can be ignored, leaving only two 

independent helicity amplitudes, Ai \ g , differing only in the gluon polarization]]. 

Because of gauge invariance, each amplitude satisfies the Ward Identity resulting from 
replacing the gluon polarization vector with the gluon four momentum. This simplifies to 
J- + p g ■ G = 0. The form factors are calculated numerically using the FF Fortran library 



[12], so the Ward Identity can be verified numerically. Rotational and Lorentz invariance 
are also checked in the same manner. 

The high-pTproduction of the Higgs boson at a hadron collider is calculated by fold- 
ing the parton-level cross section with the parton distribution functions (PDF). We use 



5 The superscript T denotes taking the transpose 
6 <7i are Pauli matrices satisfying Tr(aitjj) = 2Sij. 
7 In our numerical results, we include all the contributions. 
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CTEQ2L parton distribution functions and evaluate coupling constants at the momentum 
scale Q 2 = s. Unless otherwise stated, we use mt=175 GeV in all calculations. 

3 Numerical Results 

High-pyHiggs boson production at a 2 TeV pp collider is too small to be observed for 
all practical purposes.^ We present results only for the LHC (a 14 TeV pp collider). In 
Table I, we list the production cross section to 0(a s Gp) for several Higgs boson masses 
as well as the O(a^Gp) contribution for the transverse momentum of the Higgs boson 
pip > 30 GeV. There are separate columns for the associated production of the Higgs 
boson with 6-quarks and u, d, s, c-quarks. For all of these results, we include both the 
quark and antiquark contribution. For the associated production of H with a 6-quark or 
anti-quark, the total 0(a s Gp) cross section for pp > 30 GeV is as large as 10-20% of 
the O(a^Gp) for Mjj = 100 — 200 GeV. As the transverse momentum of the Higgs boson 
increases, the 0(a s G F ) contribution becomes relatively more important. For this process, 
the mt dependence is minimal. For instance, for Mjj = 110 GeV, the total 0(a s Gp) cross 
section of b/b+H is 16.9 fb and 17.5 fb for mt = 160 GeV and 190 GeV, respectively. For the 
associated production of H with light quark, the 0(a s Gp) cross section for pp > 30 GeV 
is never more than about 3% of the (D(a;?Gj?) for Mh = 100 — 200 GeV. The cross section 
for q + q ^ H + g is much smaller than for the corresponding process g + q/q — > q/q + H, 
having values (1.5,. 15) fb for Mh = (110,400) GeV, and will not be discussed further. 

We also studied the pp dependence of the cross section as a function of Mh- In Figure 2 

we show the cross section integrated above pp for b/b + H production to 0(a s Gp) and 

O(o^Gf) in the same Mh range as in Table I. The mean pp of the (D(a s Gp) process, for 

Pp in the range 50-350 GeV, is a slowly varying function of Mh below Mh = 180 GeV, with 

a value of approximately 120 GeV. Above Mh = 180 GeV, the mean pp ranges takes on 

the values (139,152) GeV for M H = (200,400) GeV. The mean pp* of the 0(a 3 s G F ) process 

is strongly dependent on the lower transverse momentum cutoff needed to regulate the 

Pp — ► divergence associated with the gluon progagator, and is somewhat smaller than 

that for the 0(a s Gp) process. 

Finally, we address the issue of the 0(a s Gf) process g + b/b — * b/b + H, which is a tree 
8 The rate at 2 TeV (pp) is about two orders of magnitude smaller than that at 14 TeV (pp). 
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Higgs Production, ptf > 30 GeV, m t = 175 GeV 


Mu 
(GeV) 


a(g + b/b -> b/b + £T) (fb) 


°(g + q/q^q/q + H) (fb) 


0{a a G%) 


0{a 3 s G F ) 


0(a s G 3 ) 


0(a 3 s G F ) 


110 


17.1 


125.4 


46.5 


2.2xl0 3 


140 


15.1 


88.8 


37.0 


1.6xl0 3 


180 


7.4 


61.6 


25.2 


1.2xl0 3 


400 


1.9 


29.6 


0.5 


0.7xl0 3 



Table 1: Cross section for Higgs boson production at high-p^. 

level process. To estimate the size of this cross section, we used the processes q+q — ► b+b+H 
and g + g — > 6 + 6 + i^in Pythia 5.7 |l^] with > 50 GeV. We obtained the values 
(4.5,3.0,2.0, .3) fb for M H = (110,140,180,400) GeV. In the limit that m b -> 0, there is 
no interference between the tree level process and the higher order amplitudes 9 , so any 
observed cross section is primarily the 0(c>! 3 Gf) and 0(a s G F ) processes. 

In summary, we find that to accurately predict the cross section and test the properties 
of an intermediate mass Higgs boson produced at high-py in association with b quarks 
and anti-quarks, the 0(a s G F ) contributions should be included with the 0(a 3 G^) con- 
tributions. Because of the large top quark masss, the 0(a s Gp) contributions are larger 
than the tree-level contributions of 0{a s Gp) for large . Although the production rate of 
0{a s Gp) for the associated production of the Higgs boson with a light quark or anti-quark 
is not negligible, it is only a few percent of the 0(a 3 Gi?) rate, and therefore is probably not 
distinguishable from the uncertainty in the PDF. 

4 The Electroweak Chiral Lagrangian and Non- Standard Model 
Couplings 

The process under consideration is sensitive to electroweak symmetry breaking in three 
different sets of couplings. First, there is the t-b-W vertex. The additional non-standard 
couplings can be deduced from the chiral Lagrangian [14]: 

CCh 



C = -V2Kf J t L ^b L T,+ -V2K C L c h Ll HL^ 



9 We note that the tree level amplitude and the higher order amplitudes have different helicity structure 
in the mj — > limit. 
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Figure 2: The cross section integrated above for b/b + H production to 0(a s Gp) and 
0(oIGf) for various Mh- 
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-V2K C R c t Rl »bRV+ - y/2K% c hnrft R T.- , (2) 

where £^ = -^(S^ =F for £^ = — |Tr(<7 a Xt-D /t £), and the action of the covariant 
derivative is D^T, = d^T, — igW^^T, + ig £_B^?p The matrix field £ = exp (i ^ ° ) , where 
cr a , a = 1,2,3, are the Pauli matrices, and </> a 's are the Goldstone bosons. Second, there is 
the Yukawa coupling between the Higgs boson and top quark. The most general Yukawa 
coupling of the fermion doublet F in the chiral Lagrangian is: 

C = --HFMF, (3) 

v 

where M is a 2 x 2 mass matrix. In the Standard Model, cq = 1. Third, the coupling of the 
Higgs boson and the electroweak Goldstone bosons comes from the Lagrangian: 

c = ^Hd»H -\m 2 h h 2 -V(H) 

+(%H + |# 2 )TY (p^D^) . (4) 

In the Standard Model, c\ = C2 = 1. 

To illustrate that new physics effects may enhance the 0(a s Gp) rate but not the 
0{oP s Gf) rate, we study the effects of new physics arising from the scalar sector of the 
Lagrangian, as shown in Eq.(|j). In this case, the O(a^G^) contribution is not modified. 
As shown in Ref. [ 10 1 , some models of the symmetry breaking sector allow the coefficient 



c\ in Eq.(||) to be larger than 1. (c2 is irrelevant to the processes of interest.) For instance, 



ci = a/8/3 was discussed in Ref. [|10|]. Because new physics can simutaneously modify the 
interactions of t-b-W and W-W-H, we do not intend to give predictions for any specific 
model. For simplicity, we only study the effects of new physics due to c\ in the limit that 
the SU (2)l x U(l)y gauge couplings g and g' vanish. Table II contains some of our results. 

Although these rates do not represent the true rates of the process g + b/b — > b/b + H, 
they illustrate that the rates can vary by about a factor of 2 for a heavier Higgs boson. If 
the electroweak corrections to high-p^ Higgs production are substantially modified by new 
physics at the order 0(a s Gp), then this can be observed in b/b + H production at future 
hadron colliders. 

5 Discussion and Conclusions 



Because the top quark mass is large, of the order of v, the interaction of the top quark 
and the Goldstone bosons is strong and, therefore, can be sensitive to the electroweak 
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Higgs + b/b Rate (fb) as g, g' 
p% > 30 GeV, m t = 175 GeV 


1\/T rr 

(Uevj 


Chiral Lagrangian coefficient C\ 


-.5 


0.0 


0.5 


1.0 


1.5 


110 


1.09 


0.87 


0.78 


0.83 


1.01 


140 


1.26 


0.79 


0.87 


1.49 


2.64 


180 


0.74 


0.69 


2.00 


4.67 


8.71 


200 


0.42 


0.63 


1.76 


3.82 


6.81 


250 


0.19 


0.49 


1.32 


2.70 


4.60 


300 


0.12 


0.36 


1.00 


2.01 


3.41 


400 


0.06 


0.21 


0.59 


1.20 


2.03 



Table 2: Effects of new physics on high-p-rHiggs boson production as g,g' — > 0. In the 
Standard Model, c\ = 1. 

symmetry breaking sector. For the associated production of b/b with the Higgs boson at 
high pr, the SM electroweak corrections of 0(a s G F ), involving the large top quark mass, 
is comparable to the the QCD corrections of 0(o? s Gf)- On the other hand, the associated 
production of light quarks and anti-quarks with the Higgs boson is not significant because 
no large fermion mass is involved. 

Once the Higgs boson is discovered, it is important to test whether it is a Standard Model 
Higgs boson or some other non-standard scalar particle. The cross section of the Higgs boson 
production at large transverse momentum can be sensitive to new physics which modify 
either t-b-W, t-t-H, or W-W-H vertices. Among them, only the t-t-H vertex can modify 
the 0{o? s Gf) contributions. In contrast, all of them can modify the 0(a s G F ) contributions. 
As illustrated in Table II, it is possible that the 0{a s G F ) rate is enhanced by more than 
a factor of 2 due to new physics effects. Therefore, 0(a s Gp) contributions should also be 
included when testing SM predictions. 
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Appendix I: Loop Integration 

In calculating the helicity amplitudes, one must evaluate loop integrals of the form 

J_ f d n Q{l,Q^Q»Q»} 

ITT 2 J (Q2_ m 2 )(( g + p) 2_ m 2)...- 

For triangle diagrams, X = C, and for box diagrams, X = D. 
Triangle Diagrams 

The scalar function for triangle diagrams, showing explicitly its dependent variables, is: 
C (ml ml mlplplpl) = -L J ^ _ ^ q + ^ + ^ + ^ _ ^ , 

where the internal line masses rrii are labelled by the external lines, p\ is the momentum 
flowing between the lines with masses m\ and m 2 , p 2 between mi and 7713, and p 3 = —p\ —p2 
between m% and m\. The vector integral over is 

Cup? + Ci 2 p£. 
Similary, the tensor integral over Q^Q V is 

C21KK + C 22 p£p£ + c 23 K^ + + (W- 

Box Diagrams 

The scalar function for box diagrams, similar to the triangle diagrams, can be written as 
Do(mi, m 2 , m|, m\, p\, p 2 , p\, p\, {pi + P2) 2 , {P2 + P3) 2 ) ■ The notation is an obvious general- 
ization of that for the triangle diagrams. The vector integral over Q M is 

Dup'i + D 12 p% + D 13P %. 

The tensor integral over Q^Q U is 

D 2lP »p\ + B 22 p^p\ + D 23P ^+ 
D 2 Mp v 2 + PIP2} + D2MPI + P1P3 } + D 2 MPz + P2P3} + D 2 79^. 
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Appendix II: Form Factors 

In 't Hooft-Feynman gauge, there are 20 Feynman diagrams containing Vl^-bosons and 
=t Goldstone bosons involved in the process g + b/b — ► b/b + H at 0(a s Gp) for = 0. 
As discussed in the text, the 3 diagrams involving Z°-bosons and cp° Goldstone bosons are 
negligible. Some typical diagrams are shown in Figure 1. In this appendix, we list the 
individual contributions to the form factors [cf. Eq.(Q)] from each Feynman diagram for the 
process g(q g ) + b(qs) — ► b{q\) + H{q 2 ). All momenta are defined pointing in to the Feynman 
diagram, i.e. the outgoing quark (gi) and Higgs boson (q 2 ) four momenta have a negative 
energy component. There are 12 triangle diagrams, with terms labelled T% — J- 2 q, and 8 
box diagrams, with terms T\ — T% and G\—G%- The full form factors are T = Y^f=i T an d 
G^ = Ya=i The following expressions contain the invariant masses Sij = + qj) 2 and 
s = (^i + 52 + Q3) 2 = q 2 g - We have used the relation M\y = \gv to re-express the electroweak 
coupling constants in terms of masses and the vacuum expectation value v. The process 
involving light quarks in the initial and final state can be deduced by setting mt = 0. 
The limit g, g' — ► 0, which we take to study the electroweak chiral Lagrangian, is obtained 
by eliminating all terms with an explicit M\y dependence. All triangle diagrams contain 
progagators for one fermion and two gauge or Goldstone bosons. The box diagrams fall into 
two categories, those containing two fermion and two gauge or Goldstone boson propagators 
(denoted t-t-W-W) and those containing three fermion and one gauge or Goldstone boson 
propagators (denoted t-t-t-W). 

Triangle Diagrams 



•^"10 



T 9 



-8C l2 M^/v 3 
-2C X2 m 2 t M 2 H /v 3 



T12 



2m 2 M^(-2C 00 - C l2 )/v 3 
-2m 2 M 2 ,(-C o + C 12 )/v i 



,3 



Fl3 



-8C 12 M^/v 3 
-2C X2 m 2 t M 2 H /v 3 



Ti§ 



Tis 



2m 2 M^(C 00 - C l2 )/v 3 
-2m 2 M^(2C 00 + C l2 )/v 3 
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Tyj = 2m 2 t M^(-2C 00 -4C n )/v 3 

F 18 = 2mt(-C 00 -2C u )/v 3 

F w = 2m 2 t M$ v {-2C 00 -4C 11 )/v 3 

T 2 o = 2mf(-C 00 -2C 11 )/v 3 

Note that there are no tensor contributions from the triangle diagrams, i.e. no C 2 i 

Box Diagrams 



t-t-W-W diagrams 



T 2 = 

T 3 = 

T 4 = 

= 

Q$ = 



t-t-t-W diagrams 
T 7 



-8M^(2D 27 - D 00 m 2 t + (D 22 - D 24 + D 25 - D 2% )M 2 H 

+{D n + D 25 )s + (D 24 - D 25 )s 12 - (D n + D 12 - D 25 + D 2% )s 23 )/v 3 

-2m 2 t M 2 H {2D 27 - D 00 m 2 + (D 22 - D 24 + D 25 - D 26 )M 2 H + 

(£>13 + D 25 )s + (D 12 - D 13 + D 24 - D 25 )s 12 - (D 25 - D 2& )s 23 )/v 3 

-2m 2 M^(-2D 27 + Diss + (2D Q0 + D 12 - D 13 )s 12 )/v 3 

2m 2 M^{2D 27 + (Doo + D u )s + (-2D 00 - D u + D 12 )s 23 )/v 3 

-16M^((-D n - D 24 )q» + (-D 12 - D 22 )q^ + {-D 12 - D 26 )q%)/v 3 

-Am 2 M 2 H ({-D l2 - D 24 )q? + {-D 12 - D 22 )q^ + {-D 13 - D 26 )q^)/v 3 

-4m?M^((-2A)o - 2D U - D 12 - D 24 )q» 

+ (-2A)o - 3D 12 - D 22 )q^ + (-3D 13 - D 26 )q^)/v 3 

4m 2 M^{(-2Doo - 2D U + D 12 + D 24 )q» 

+ (-D 12 + D 22 )q^ + (-D 13 + D 26 )q%)/v 3 

4m 2 M^{Do m 2 + {-D 22 + D 24 - D 25 + D 26 )M 2 H - (2D U + D 25 )s 
+ {-D 2 4 + D 25 )s 12 + (-Doo + 2D U - 2D 12 + D 25 - D 26 )s 23 )/v 3 
2mj{Doom 2 + (-D 22 + D 24 - D 25 + D 26 )M 2 H - (2D 13 s + D 25 )s 
+(-D 00 - 2D 12 + 2D 13 - ^4 + D 25 )s 12 + (D 25 - D 26 )s 23 )/v 3 
4m 2 M^{Doom 2 + {-D 22 + D 24 - D 25 + D 2% )M 2 H 
+ {D 00 - D n + D 13 - D 25 )s + (-D 00 + D 12 - D 13 - D 24 + D 25 )s 12 
+(D U - D 12 + D 25 - D 26 )s 23 )/v 3 
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2m\{D m m\ + {-D 22 + D 24 - D 25 + D 26 )M 2 H 

+ (D 00 + D u - D 13 - D 25 )s + (-D 12 + D 13 - D 24 + D 25 )s 12 

+ (-D 00 - D u + + D 25 - D 26 )s 23 )/v 3 

8m 2 t M^{(3D u + 2D 24 )q» + (D 00 + 3D 12 + 2£> 22 )g£ 

+ (D 00 + 2D 12 + D 13 + 2D 26 )q%)/v 3 

4mt((D 00 + D u + 2D 12 + 2D 24 K + (A)0 + 3£>i 2 + 2Z? 22 )^ 

+(3D 13 + 2D 26 )^)/t; 3 

8m$M$y((D u + 2Z? 24 )< + (D12 + 2£> 22 )«£ 

+ (-^00 + 2Di 2 - L>i3 + 2D 26 )q^)/v 3 

Am A t ({-D m - D n + 2L> 12 + 2£> 24 K + (£>i 2 + 2D 22 )^ 

+(D 13 + 2D 26 )q$)/v 3 
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